Competing orders II: the doped quantum dimer model 
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We study the phases of doped spin S — 1/2 quantum antiferromagnets on the square lattice, 
as they evolve from paramagnetic Mott insulators with valence bond solid (VBS) order at doping 
5 = 0, to superconductors at moderate 5. The interplay between density wave/VBS order and 
superconductivity is efficiently described by the quantum dimer model, which acts as an effective 
theory for the total spin S = sector. We extend the dimer model to include fermionic 5 = 1/2 
excitations, and show that its mean-field, static gauge field saddle points have projective symmetries 
(PSGs) similar to those of 'slave' particle U(l) and SU(2) gauge theories. We account for the non- 
perturbative effects of gauge fluctuations by a duality mapping of the S = dimer model. The 
dual theory of vortices has a PSG identical to that found in a previous paper (L. Balents et at, 
|cond-mat/0408329 ) by a duality analysis of bosons on the square lattice. The previous theory 
therefore also describes fluctuations across superconducting, supersolid and Mott insulating phases 
of the present electronic model. Finally, with the aim of describing neutron scattering experiments, 
we present a phenomenological model for collective 5=1 excitations and their coupling to superflow 
and density wave fluctuations. 



I. INTRODUCTION 

In a previous paper^ (hereafter referred to as I), we de- 
scribed the physics of two-dimensional superfluids in the 
vicinity of a localization transition into a commensurate 
Mott insulator. Our primary focus was on bosons on the 
square lattice, and we demonstrated the role played by 
the projective transformations (the PSG) of the vortices 
in the superfluid under operations of the square lattice 
space group. The present paper will extend our analysis 
to paired electron models on the square lattice, with the 
aim of justifying applicability to the cuprate supercon- 
ductors. A simple physical argument based on universal- 
ity was given in I clearly validating that analysis when 
the superconducting state neighboring the Mott insula- 
tor is an s-wave state with a large quasiparticle gap, and 
it was claimed that the same description also holds for a 
similarly gapped "strong pairing" d-wave state. In this 
paper, we back up this claim with detailed microscopic 
analysis of a set of models that incorporate local singlet 
formation. Specifically, we find that most of the results 
of I can be applied essentially unchanged, with the PSG 
determined by the particle density of Cooper pairs. We 
will also discuss the extent to which our results can be 
applied to "weak pairing" d-wave superconductors with 
gapless nodal quasiparticles. 

Most of our analysis here will be carried out in the 
context of the quantum dimer model. This model was 
proposed by Rokhsar, Kivelson, and FradkinSii as an ef- 
fective theory of the S = excitations of the insulat- 
ing, paramagnetic phases of square lattice antiferromag- 
nets. For generic parameters, the insulating dimer model 
has a ground state which breaks the symmetry of the 
square lattice space group with the development of va- 
lence bond solid (VBS) orderiSii^. The dimer model was 
also extendedSii to allow for S = Q charged excitations: 



the Hilbert space was expanded to include bosonic holes 
with density 8, and it was argued that the ground state is 
a superfluid at finite 5 (the pairing symmetry of this su- 
perconducting state is somewhat unclear in these works, 
and will be clarified below). The doped quantum dimer 
model is therefore ideally suited to our primary purpose 
here of describing the evolution from the VBS Mott in- 
sulator to the superfluid as a function of increasing 5. 

We will begin in Section ^ by reviewing the 
formulation^'^ of the quantum dimer model as a compact 
U(l) gauge theory (with gauge field A) in the presence 
of a background static matter field of charges ±1 on the 
two sublattices of the square lattice. Upon doping, the 
holes appear as two S = Q dynamic bose fields hie {i is 
a square lattice site index and the index ct = ±1 labels 
the sublattice upon which the boson primarily resides). 
These bosons have charges a under the compact U(l) 
gauge field A. 

In Section lllll we show how the above theory of the 
spin 5 = sector can be extended to include 5=1/2 
fermionic excitations fis (s =tii is a spin label). These 
fermionic 'spinous' also carry unit staggered charges 
under the compact U(l) gauge field A. While these 
fermionic excitations are gapped and confined in the un- 
doped insulator, the structure of our theory allows these 
fermions to become gapless quasiparticles at nodal points 
in a possible d-wave superconducting ground state at fi- 
nite 5: they then become the conventional 5 = 1/2 Bo- 
goliubov quasiparticles of this superconductor. With this 
construction, we can directly discuss the pairing symme- 
try of the superfluid states of the doped dimer model. We 
will also consider collective 5=1 'triplon' excitations in 
Section El 

So far we have mainly defined the degrees of freedom of 
our theory on the lattice scale, but this generally does not 
specify its ultimate ground state or the quantum num- 
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bers of its quasiparticles. In particular, we are consider- 
ing a strongly interacting 2+1 dimensional compact U(l) 
gauge theory, and the non-perturbative fluctuations of 
A must be treated with great care. Generally, for the 
dimer models we consider here, these fluctuations drive 
confinement at long scales. Nevertheless, it is useful to 
initially neglect gauge fluctuations and to consider mean- 
field saddle points with a fixed background A. Such sad- 
dle points describe (unstable) 'spin liquid' states, which 
may be considered the progenitors of the superconduct- 
ing and/or charge ordered Mott phases which appear as 
ground states. A subsequent treatment of gauge fluc- 
tuations describes this evolutionary descent via confine- 
ment or the Higgs mechanism. Remarkably, we will find 
in Section Fill CI that, at least at the level of PSGs, the 
spin liquid saddle points of the dimer model are related 
to particular (gener alized ) saddle points of earlier slave- 
particle U(l) (Refs. or SU(2) (Ref. |n|) gauge the- 
ories. The particular saddle points appropriate to dimer 
models, however, are distinct from those advocated as 
candidate spin liquid states in the slave particle work. 
They have appeared - without much attention - in Wen's 
PSG classification of symmetric spin liquidsi^ (the spe- 
cific identifications of these spin liquids appear at the 
end of Section IIII ("ll . Our work thus resolves one ques- 
tion posed in Rcf. 16, as to how the (putative) spin liquid 
states of the quantum dimer model fit into this classifica- 
tion scheme. Beyond classification, the physical interpre- 
tation we shall give to our dimer model saddle points, and 
our subsequent non-perturbative treatment of gauge fluc- 
tuations, lead us eventually to conclusions on the physical 
properties of our theory which differ from those reached 
in Refs. lnllTl 

There is one situation in which it is clear that the A 
gauge fluctuations are relatively innocuous. This is the 
moderate S state in which the bia- are strongly condensed 
and the square lattice symmetry is fully preserved: the 
condensate acts like a Higgs field which quenches A fluc- 
tuations. Such a state is a superconductor, and we will 
show in Section Fill Dl that microscopic energetic consid- 
erations and the PSG constrain the pairing symmetry of 
the superconductor; d- and s- wave (and other) pairings 
emerge under suitable conditions. 

So what is the fate of this putative d-wave supercon- 
ductor as S is reduced? As noted above, we know from 
Refs. [in that strong A fluctuations produce a VBS in- 
sulator at (5 = 0, and so a careful treatment of A fluc- 
tuations must surely become increasingly important as 5 
is reduced. In previous workiSii^ii^ii^ on U(l) slave par- 
ticle and related theories, the evolution between a flnite 
d d-wave superconductor and a VBS insulator at S — 
was addressed in a mean-field theory (similar to that just 
discussed above) but which also allowed for breaking of 
the square lattice symmetry by the development of den- 
sity wave/ VBS order coexisting with superfluidity. This 
had the advantage of explicitly displaying the evolution of 
the gapless fermionic nodal quasiparticle excitations with 
decreasing S: the nodes present at larger S disappeared 



(due to collisions in pairs at the (reduced) Brillouin zone 
boundaries) at a critical doping S — Sf, larger than the 
doping Sc at which the Mott transition occurs. For small 
enough S < dc, a fully-gapped insulating VBS state was 
obtained. However, this work treated gauge or VBS fluc- 
tuations in a cavalier manner, and the VBS order was 
static and likely overestimated in magnitude. 

In the application to the cuprates, we have in mind 
that these insulating or supersolid states obtained in our 
theory at finite doping are models of the situation near 
doping S = 1/8. There is much experimental evidence 
{e.g. in Lai.6-a;Ndo.4Sra;Cu04 in Ref. IT^ that exactly 
at (5 = 1/8 the superconductivity is strongly suppressed, 
and perfect long-range order with a period of 4 lattice 
spacings develops (as is also the case in Fig^. 

In Section IIVI we will address the 5 driven evo- 
lution from a VBS insulator to a superconductor 
in an approach which fully accounts for the non- 
perturbative effects of the fluctuations of A. No previous 
analysisi2iii*i2iiii^iiSii& has accounted for strong gauge 
fluctuations in this regime. However, we will only be able 
to do this in the context of the S' = model of Sect ion ITTI 
non-zero S excitations will be neglected in Section IIVI 
While this is probably safe within the phases obtained 
therein, neglect of gapless nodal S' = 1/2 fermionic exci- 
tations is likely invalid at the quantum critical points. A 
full treatment which synthesizes the results of Sections llIII 
and IIVI and thus accounts for both strong A fluctuations 
and gapless 5' = 1/2 fermions at quantum critical points 
is not available, and remains a challenging and important 
open problem. 

The analysis of the 5 = sector in Section llVl is carried 
out by a duality mapping into a theory of vortices. The 
resulting effective theory for the vortices has a structure 
similar to the models considered in I, and can therefore 
be similarly classified by an analysis of the action of the 
PSG on the vortices. With the knowledge of the PSG, 
many of the results of I can be directly applied here. We 
will also present in Section llVBI a mean-field analysis of 
the dual vortex theory which results in a phase diagram 
of the dimer model as a function of S. It is important to 
note that this dual mean-field theory accounts for non- 
perturbative confining effects of the A fluctuations: these 
are built into the duality mapping. 



II. QUANTUM DIMER MODEL 

This section describes the doped quantum dimer 
modelSi, viewed as an effective theory of the spin S = 
sector of doped quantum antiferromagnets. The resulting 
theory will be a compact U(l) gauge theorjtL^ coupled to 
two classes of matter fields: (?) static charges of ±1 on 
the two sublattices, and {ii) dynamic bosons of charges 
±1 with a total physical density S (note: bosons of both 
gauge charges contribute a positive number to the phys- 
ical density). 

The Hamiltonian of the quantum dimer model can be 
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FIG. 1: Elementary moves of the quantum dimer model, (a) 
The resonance of two dimers around a plaquette, which occurs 
with amplitude J in Eq. I|2.3|l . (6) Hopping of a hole with the 
amplitude t term in Eq. 12.71 . The hole (represented by the 
open circle) prefers to reside on one sublattice. This hole 
is on the "wrong" sublattice (with high energy cost) in the 
intermediate state above, which only acts as a technical tool 
for effectively generating the same sublattice hopping in the 
low energy limit. The constraint in Eq. (12.811 is obeyed on all 
sites in all the states shown. 

written as 

Hqd = Hd + Hh + Ht (2.1) 

The first term, Tid involves a diagonal energy and a res- 
onance of dimers on a plaquette, as illustrated in Fig^. 
Following Ref. 0, and the notation of Ref. H, the dimers 
are represented by the conjugate phase Aia and integer 
number Eia variables which reside on the links of the 
lattice: 

[Aia,Ejp\ = iSijSap ■ (2.2) 

The dimer number on a bond is rjiEia oriented to the 
right or the top {rji = +1 or -1 on the two sublattices). 
The Hamiltonian for dimers on a plaquette is 

'^d^jJ2^^c, + jJ2 (e«/3Aa A/3) (2.3) 

ia i 

Note that the co-efficient of the cosine has the "wrong" 
sign. This is a consequence of the microscopies of the 
electrons in the t-J model, if one adopts the natural sign 
convention that a singlet dimer is created by an electron 
pair-field operator (see Eq. 1)3. and a given dimer cov- 
ering corresponds to a product of such operators acting 
on the vacuum. 

The next two terms, Tih and Tit, describe the hop- 
ping of holes in the doped dimer model, as illustrated in 
Fig ^p. To maintain the description in terms of near- 
est neighbor dimers, the holes must hop between sites 
on the same sublattice. Consequently, the model has 
two species of holes, one for each sublattice. This same- 
sublattice hole hopping model is somewhat inconvenient 
to work with under the duality mapping, and so we ex- 
pand the Hilbert space of the model by additional high 
energy states to allow us to express all hole motion in 



terms of nearest neighbor hopping: we allow each hole 
species to temporarily hop onto the "wrong" sublattice 
in a state which has a large energy. This hole will quickly 
hop back onto the "right" sublattice, and the two hop 
process will have achieved the same sublattice hopping 
we are seeking to model. This process is illustrated in 
FigCb; the same strategy was used by Moessner et alm& 
in their study of the doped dimer model. We represent 
the two species of holes by the rotor phases 0i- and 
the conjugate number variables n^-. The commu- 
tation relations are 

[(t>icr,nja'] = iSijScrcr' (2.4) 

where a = +, —■ These holes reside primarily on the -I- 
or — sublattices and will have opposite gauge charges. 
The potential energy of the holes is described by Hh'- 

H/t = ^ {h^o■ - H - arjiWf + ^ Ki^fHaflja ■ 

(2.5) 

Here ?7 > is the hole self repulsion, 2H is the aver- 
age hole density preferred by the chemical potential (we 
are working at fixed chemical potential, and the actual 
density, 5, in the ground state of the Hamiltonian may 
be different), rji = ±1 on the two sublattices, and 2UW 
is the energy penalty for a hole to be on the "wrong" 
sublattice. For W > 0, the -I- holes are preferred on the 
rji — 1 sublattice, and the — holes on the other. The off- 
site interaction energy Aij is assumed to be repulsive, and 
is required to stabilize insulating states away from zero 
doping. We will discuss the case of general , although 
our numerical results later will be restricted to short- 
range Aij] long-range Coulomb interactions lead only to 
minor modifications. The total density of holes on site i 
is given by 

S^ = J2(n^a) . (2.6) 

The nearest neighbor hopping Hamiltonian for the two 
hole species is 

Ht^-tY^ cos (a„(^.,^ + crAa) ■ (2.7) 

iaa 

Finally, there is a local constraint, which commutes 
with all the terms in the Hamiltonian, 

^aEia + ^ Cn%a- = Tji ■ (2.8) 
a 

This constraint ensures that there is either a hole on each 
lattice site or a dimer emerging from it. ft also allows 
for a configuration with two dimers emerging from that 
site if the extra gauge charge is compensated by an ad- 
ditional hole from the "wrong" sublattice. This is the 
intermediate state in Fig which compared with the 
other two states in that figure requires the excitation en- 
ergy V/2 ^ UW. We see that for sufficiently large V 
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we can safely set W equal to zero and still pay a price 
for having a hole on the "wrong" sublattice. Choosing 
both V and U sufficiently large also excludes unphysical 
states with even larger excitation energies. Note that all 
the states in Fig obey Eq. Ij'i.Sfl . 

We have already noted that the dimer resonance term 
in Eq. (|2.3|l has the opposite sign from that usually found 
for the Maxwell term in lattice gauge theories. We will 
see below that this can have significant consequences for 
the physics, and hence the structure of the dual theory. 
Physically, the kinetic term in Eq. H2.3|l favors an average 
TT flux per plaquette in the gauge field Ata- The hole 
propagation term, Eq. 1)2.711 . however, clearly allows holes 
to achieve lower kinetic energy if this flux is expelled. 
Since the gauge field is dynamical, the system will choose 
for itself what flux penetrates each plaquette on average, 
based on a balance of these two competing energies, and 
as pointed out by Kivelson and Fradkin'', the flux will 
penetrate at 5 = and be expelled for large enough 6. In 
a strongly confining phase, as occurs for 6 = 0, however, 
this fiux has no physical significance, since the gauge field 
is strongly fiuctuating. 

Mathematically, one may choose to view either the zero 
or TT flux states as vacua, provided fluctuations of this 
fiux are included. In particular, one may attempt to fix 
the sign "problem" (J > 0) by shifting the gauge field 
Aia on certain links by tt. However, this has the effect 
of inducing a tt flux per plaquette in the motion of the 
(/fio- holes in Eq. 1)2.711 . If we were to perform such a 
transformation, Tid and Tit would be replaced by 

Wd = |E^»-'^EcoK^"/3^"-^^'3)' (2-9) 

i i 

H't = -tY^ cos (^Aa$ia + aAia + 27r.g,„) , (2.10) 

iaa 

and the holes now move in a magnetic field gi^ which 
obeys 

e^^AA^g^A = hS^r (2-11) 

with h=l/2. 

For (most of) our purposes, it will generally not be 
advantageous to make such a shift. This is physically 
because our description of the Mott transition approaches 
the critical point from the superconducting side, in terms 
of the condensation of vortices of the superconducting 
state. As discussed above, in the large i5 region of the 
dimer model where superconductivity occurs, zero flux is 
favored. Hence this is a more natural point from which 
to embark upon our (dual) investigation. 

In the following section, however, in which we discuss 
quasiparticle physics and pairing symmetry, we will con- 
sider the physical content of this 7r-fiux background. Our 
intuition to expand around the zero flux configuration is 
there confirmed, as we show that the 7r-flux background 
is inconsistent with d-wave (or s-wave) pairing in the su- 
perconductor, and applies instead to a yet more uncon- 
ventional state. 



III. ADDING FERMIONS AND SPIN TO THE 
DIMER MODEL 

In the context of the cuprates, a shortcoming of the 
present approach - intrinsic to the dimer model - is the 
neglect of fermionic spin-carrying quasiparticles. To un- 
derstand the nature of this difficulty more deeply, we 
discuss the results of a natural extension of the dimer 
approach to include these = 1/2 fermionic excitations. 
The dimer model can also be extended to include collec- 
tive S = 1 'triplon' excitations, but we defer discussion 
of these to Section EI until after we have performed a 
duality analysis of the S = dimer model in Section Hvl 

Readers mainly interested in the duality analysis of 
the S = dimer model at finite doping may skip ahead 
to Section llVI and return to the present section upon a 
second reading. 

The strategy will be to first understand in more detail 
the nature of the insulating states of this model when 
undoped. Because in this model at half-filling all the 
"matter" (holes and fermions) is gapped, the Polyakov 
argument implies that the U(l) gauge theory must be in 
a confining phase. Nevertheless, it is instructive to con- 
sider formally the limit of weak gauge fluctuations, which 
in practice applies for large J - of either sign. Then one 
may view the insulator as a gapped U(l) state plus gauge 
fluctuations, the latter - more specifically monopole pro- 
liferation - ultimately leading to confinement on long 
length scales, and with it, density wave order. This view 
has the advantage that the superconducting state can be 
seen as a "Higgs" phase obtained from the U(l) state. 
Further, the insulator can be characterized by a PSG, 
defining the manner in which the physical symmetries 
are realized in the gauge theory. Here we will consider 
two such PSGs, which correspond naturally to the phases 
obtained for large positive and negative J. 

We can better understand in this manner the nature of 
the superconducting and insulating states in this model. 
We will see that the former should be understood as a 
"molecular" or "strong pairing" d-wave phase. The latter 
can be seen to arise from either of two natural underly- 
ing U(l) spin liquid states. These two U(l) states are 
distinct, and characterized by different PSGs. They are, 
however, both unstable at low energies to the various 
charge ordered states discussed herein. 

We proceed with a multi-pronged attack. First, we de- 
scribe the symmetry of the bosonic dimer -|-hole Hamil- 
tonian of Sec. ^ We consider the two distinct cases in 
which the ground state has an average zero or tt gauge 
fiux per plaquette. These two cases give rise to distinct 
PSGs which describe different putative insulating states. 
With this in hand, we discuss the nature of the pairing 
symmetry of superconducting states obtained from this 
model. Next, we will augment the dimer-|-hole Hilbert 
space by allowing each site to contain an unpaired elec- 
tron (or spin), in addition to the end of a dimer or the 
hole allowed up to now. Physical arguments give the gen- 
eral form of the Hamiltonian with such unpaired electrons 
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included. This form is further restricted by combining 
the PSGs of the dimer+hole model and the requirement 
that physical electrons transform properly. We discuss 
how this analysis is related to slave particle formulations 
of Mott insulating states and Wen's analysis^^ of their 
PSGs. With the full form of the dimer+hole+electron 
Hamiltonian determined, we confirm the symmetry anal- 
ysis of the superconducting states by a direct investiga- 
tion of the effective Bogoliubov-deGennes quasiparticle 
Hamiltonian. 



A. Projectively realized symmetries of the doped 
dimer model 

It is clear that, even without including fermionic exci- 
tations, the dimer model in Sec. ^already requires a non- 
trivial projective realization of certain symmetry opera- 
tions of the space group. This is so because the staggered 
static gauge charges (and the staggered potential term W 
in Eq. (|2.5I) ) naively break lattice symmetries. To over- 
come this, unit translations and 7r/2 rotations about a 
dual lattice site must be accompanied by the transfor- 
mations n+ h-, 0+ <-!■ 0_, A —^ —A, E — > —E to 
leave the dimer model invariant. Actually, it will prove 
useful to discuss projective symmetries with respect to a 
fixed average background U(l) flux. This will make con- 
tact with Wen's notion of projective symmetry analysis 
of slave particle saddle points. We will consider the two 
cases of zero and tt average flux per plaquctte. 

For compactness, we will in this section denote the 
boson operators by 

6,, = e-*^-. (3.1) 

Further, it is sometimes helpful to take the limit W oo, 
corresponding to a "microscopic" dimer-|-hole model in 
which two holes of course cannot occupy the same site. 
In this case the type 1 holes live only on the A sublattice, 
and conversely the type 2 holes live only on the opposite 
B sublattice. One then defines a single boson field via 

^ _ I 6"'"^'+ for i G A sublattice ,^ 
' ~ I e-"-^^- for i e B sublattice ' 

i.e. hi = bia if i is on the a sublattice. For W — > oo, 
all finite energy states can be described in terms of the 
bi operator. For < oo, this is not the case, but bi 
defined as above is still useful and represents properly 
the low energy hole degrees of freedom. 

1. Zero flux background 

First, and simplest, is the case of zero flux, the naive 
choice for J < (note that which saddle point is chosen 
by the system is however a dynamical question, which 
will in general depend upon all the degrees of freedom. 



and may vary as other parameters in the Hamiltonian are 
varied). In this case, one can easily see that the only non- 
trivial required aspect of transformations of the bosons is 
the interchange of the two boson flavors - and associated 
change of the sign of gauge fields - for space group oper- 
ations that interchange the two sublattices. Specifically, 
these are unit translations (or generally translations by 
an odd number of total lattice units), tt/2 rotations about 
a dual lattice site, and inversions through a row/column 
of the dual lattice. However, there is a mathematical 
freedom to compose each space group operation with a 
global U(l) phase rotation. That is, for an operation 
which preserves the two sublattices, one may take a pri- 
ori 

^..^^'GW.e'^^''^- (3.3) 

Similarly, for an operation that interchanges the two sub- 
lattices, 

^ foG(^),-.e'>^(«) . (3.4) 

Using the definition Eq. (|3.2|l . both these equations are 
replaced simply by 

h ^ bGi^)e'^'•''y (3.5) 

Here G{i) is the site i is mapped to under the space 
group operation G, and x(G) is the phase, which we take 
as the same for both boson flavors, and can be chosen 
differently for each G. Note that this phase rotation is 
not a U(l) gauge rotation (which would be opposite for 
the two boson flavors) unless e*-^*^'^-' — ±1. Since the 
hole fields are not themselves physical, it is not a priori 
obvious which value of x(G') is the correct one. 

An important constraint on x{G) is that any sequence 
of space group operations which compose to the triv- 
ial identity operation in the original space group should 
compose to a pure gauge U(l) transformation in the 
PSG. The square of an inversion or the fourth power of 
a tt/2 rotation must be a pure U(l) gauge transforma- 
tion. Hence the phase e''*^'^^-' = 1 if / is an inversion, 
while e^'^*^^^ = lifi?isa7r/2 rotation. Other similar 
constraints can be applied with more detailed considera- 
tions. 

For the doped dimer model, we can further refine this 
by using the physical interpretation of the hole opera- 
tors. Consider the action of the physical singlet pair field 
operator on a bond: 

= 4e,,,cl+,^_^„ (3.6) 

with a — x,y. Clearly, this operator, when acting on 
a state in which the bond in consideration is occupied 
by two holes, annihilates the holes and creates a dimer. 
Hence we require 

=/iafozfoz+e„e"'-^"". (3.7) 

Here we have introduced two unknown coefficients, fj,x , fJ^y 
(these are assumed independent of i, as can be argued is 
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generally true). To determine them, we use the known 
transformation properties of the pair field defined from 
Eq. First, it is invariant under inversion through 

the row/column of the dual lattice dividing the two sites 
of the pair. This implies e^*^^^^ = 1 (actually this is true 
for all inversions). Then e^^^-^^ = ±1, and for such values 
the phase rotation in Eq. H3.4|) is a gauge U(l) rotation. 
Thus we can always choose 



(3.8) 



Second, consider a 7r/2 rotation. Under rotation around 
site i, Eq. (|3.6|l implies "ifix '^iy, '^iy "^fi-g^^x, while 
from Eq. H3.7|) . one finds 



My 



Hence 



My 



±1. 



(3.9) 
(3.10) 

(3.11) 



No further requirements on x(^) and fia follow from 
symmetry alone. The two choices in Eq. (|3.11|l simply 
represent different possible PSGs. Note that this choice 
is directly related to the pairing symmetry of the super- 
conducting state. If fj,x = IJ'y, then a uniform condensate 
(bi) = constant represent an s-wave paired state, while 
for fj,x = —fJ-y, the same condensate represents dx2_y2- 
wave pairing. More generally, using Eq. (|3.7|) . one can 
show 



'^iy'^i+e^.y 
ix'^ i+ey ,x 



(3.12) 



where the curl in the exponential on the right hand side 
is the gauge flux through the plaquette with site i at 
the lower-left corner. Hence for the zero background flux 
states we consider in this subsubsection, using Eq. H3.11|l . 
one sees that 



^iy^i+e^,y 



1. 



(3.13) 



'^ix^i+ey,x 

This allows for both s and d-wave pairing. 

2. TT flux background 

Let us now turn to the more complicated case in which 
half a quantum of background gauge flux pierces each 
plaquette of the direct lattice. We write 



A;, 



A, 



A,. 



(3.14) 



where Aia is a classical c-number background vector po- 
tential representing tt flux per plaquette. An extremely 
convenient and symmetrical choice for this section is 



•A-i, 



Aiy 



-111 



(3.15) 



One can readily see that, with this choice, the Hamil- 
tonian remains invariant under the same translational 
operations as for zero background flux, i.e. for a unit 
translation. 



: 



(3.16) 



and Aifj — > ~Ai+e^,f3 etc. Rotations about dual lat- 
tice sites, and inversions are also unchanged. Rotations 
around direct lattice sites are, however, non-trivial. For 
a 7r/2 rotation about a direct lattice site, one requires 



hi 
bi2 



(3.17) 
(3.18) 



where 



1 for (i„ij,) = (0,0)(mod2), 

-i for (i^,ij,) = (1,0) (mod2), 

-1 for iix,iy) = (l,l)(mod2), 

i for (ij;,?^) = (0, 1) (mod2). 



(3.19) 



Note that, since Q is real on the A sublattice, one has 
simply 



CbR(^)■ 



(3.20) 



We can now repeat the symmetry analysis done above 
for zero background flux. As before, assuming the form of 
Eq. (|3.7|) , since the transformation law for inversions has 
not changed, one still obtains Eq. (|3.8|) . Under rotations, 
however, one obtains the conditions 



My 



(3.21) 



Note the important factor of i in Eq. H3.21|l compared 
to Eq. (|3.11|l . This implies very different pairing states 
if the bosons condensed in the tt flux background. For 
instance, if the holons condense uniformly, (bi) = const., 
then '^iy — T'i'^ix — const.. This is incompatible with 
either d- or s-wave pairing, and indeed with rotational 
invariance (in a d-wave superconductor, rotational sym- 
metry is unbroken for gauge-invariant observables). 

One way to see this is to consider the ratio in Eq H3.12|l . 
which becomes in this case 



'^iy'^i+e^,y 



(3.22) 



irrespective of the spatial pattern of boson condensation. 
In more conventional terms, the superconducting states 
with holons condensed in the 7r-flux background have a 
pair wavefunction of (1 ± i)sx2+y2 -I- (1 =p i)dx2-y2 form. 
Again it is clear that such a complex linear combination 
of extended s- and d-wave pairing breaks fourfold rota- 
tional symmetry. 
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B. Fermions 

Having understood the possible PSGs of the 
dimer+boson model, we are prepared to discuss fermionic 
excitations. In principle, there are many distinct ways in 
which the fermions might transform consistent with the 
above PSGs. There is, however, clearly a physically pre- 
ferred and simple choice. We proceed by allowing a site 
to be occupied not only by a hole and an end of a dimer, 
but also by an unpaired electron. This enlargement of 
the Hilbert space is represented by a spin-1/2 fermion 
mode on each site, described by canonical fermion cre- 
ation/annihilation operators fig, fig (s is a spin-1/2 in- 
dex). The original model is to be viewed as the fermion 
vacuum. Clearly, the electron annihilation operator may 
act on a site with such an unpaired electron to create a 
hole and remove the fermion. Hence we require 

c^s ~ + • • • (3.23) 

the ellipses are included to emphasize that of course mi- 
croscopically the electron annihilation operator can cre- 
ate other states. For instance, it may act upon a site 
which is part of a dimer to annihilate the dimer, create a 
hole on the site acted upon and create a fermion on the 
other site shared by the dimer; such terms would be of 
the form 

c,s bU-^'"ess'fl,^^,, + -- - ■ (3.24) 

Such additional terms are of quantitative importance in 
analyzing a microscopic dimcr+holc+fermion model. In 
any case, all such terms are such that the constraint in 
Eq. H2.8|l is now expanded to 

Aa^^a+^firi,, -|-r;,^/2,/,, =77^. (3.25) 

(7 S 

This constraint makes it clear that the fig fermions carry 
staggered U(l) gauge charge. 

The single term in Eq. (|3.23|) . however, suffices to de- 
termine the transformation properties of the fis oper- 
ators, from the requirement that the physical electron 
operator Cis is a scalar. First, we see that the fis fermion 
carries a staggered gauge charge of opposite sign on the 
two sublattices. Next, we note that with zero background 
flux, bi is a scalar apart from the x(G') factor discussed 
in the previous section. Hence 

fis e*-^'-*^-'/(3(j)5, for zero background flux. 

(3.26) 

With the TT flux background, Eq. H3.26|l holds for trans- 
formations other than the ■n/2 rotation about a direct 
lattice site. For this, one has 

^ e'^(^)C*/i?Ws, for ^ background flux. (3.27) 

At this point we note that the similarity between 
Eq. H3.23|l and the means of introducing slave particle 
"holon" and "spinon" operators in the U(l) mean field 




(a) 



f ® ^ ® f 

(b) 




FIG. 2: Extension of Fig to include 5 = 1/2 excitations 
in the dimer model, (a) Two 5* = 1/2 fermions forming a 
singlet valence bond, representing Eq. (13.2811 . (6) A holon 
and a, S = 1/2 fermion exchanging positions, representing the 
first term in Eq. 13.301 . (c) Motion of a S = 1/2 fermion 
involving rearrangement of valence bonds, representing the 
second term in Eq. 13.301 . 

theory of the t — J model is not coincidental. Indeed, the 
universal aspects of this discussion can all be recovered 
as appropriate saddle points of that approach. They can 
also be recovered from the SU(2) mean field theory, the 
two approaches being interchangeable in this instance, as 
will be returned to below. 

Having discussed the transformation properties of the 
fermions, we proceed to construct their Hamiltonian - 
assuming initially zero background field. Physically, the 
simplest allowed process is one in which a dimer converts 
to a singlet fermion pair (see Fig 01: 

i 

+ Y.5^^"'^^-fles,fl,^^,+Ii.c.. 

Here v represents the energy cost to introduce a single 
fermion, and the sums over the repeated spin indices s, s' 
are implicit. Taking w — > oo recovers the original dimer 
model. The coefficients 5x,5y appear as "pair fields" for 
the fis fermions. They are constrained by symmetry to 
satisfy 

5y = e~'^^^''^'^ 5x = ±^2;, for zero background flux. 

(3.29) 

We would also like to allow for hopping processes for 
the fermions. A fermion can hop between neighboring 
sites by exchanging places with a hole, or a fermion can 
hop between second neighbor sites if an intermediate 
neighboring site is connected to the destination site by a 
dimer, the dimer reconnecting after the hop to the origin 
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site. These processes are represented by 



(3.30) 



+ f 



isJ i+w,s^ i,i+'u^' H.C. 



where 



and w = ix + eyjW — ix — Cy. The complicated- looking T 
factors are simply the uniform superposition of the gauge 
connections (exponential of the discrete line integral of 
the gauge field) taken along the two shortest paths con- 
necting diagonal sites. The equal superposition of both 
paths is required to preserve inversion symmetry through 
the (11) and (11) axes of the square lattice. We give these 
t' terms primarily for completeness, as they will not have 
an important role in much of what follows. 

By examination of Eqs. (|3.28|1 . (|3.30l) . one sees that 
the fis fermions hop unassisted (i.e. without accompa- 
nying holes) only between sites of the same sublattice, 
and experience pair fields on links connecting sites of op- 
posite sublattices. This is indeed required by the stag- 
gered gauge charge of the fis variables. It is thus some- 
times convenient to define new variables via a particle- 
hole transformation: 



fis for ?7,j = +1, 
^ss'fL for = -1 



(3.32) 



These dis variables have a constant unit gauge charge on 
every site and transform as global SU(2) spinors. They 
should thus be identified as "spinons" . In these variables, 
one finds that the constraint in Eq. H3.25|l now takes the 
form 

A„^,„ + 6,1+ 6,+ - blb^_ + J2 dU^s = 1, (3.33) 



while the Hamiltonian Ti./ = T^p + TL^ has the form 

i&A,a 

+ E ^^+iAs^ss4+,^,s' + H.c] , (3.34) 

and we have dropped the t' terms for brevity. Here the 
pairing terms for the fis fermions have become hopping 
amplitudes: 

_ for m = +1, 

~ ^ -5*^ for ?7, = -1 . ^'^■'^^> 



From Eq. (|3.34() one can clearly see that, if no holes 
are present (or if the bi excitations are gapped), then 
the t term can be neglected, and the spinons dis obey 
a simple tight-binding model. They indeed carry unit 
gauge charge, since they are minimally coupled to Aia- 
We expect that the confining effects of A fluctuations 
will produce VBS order in the insulator, and also bind 
the spinons into pairs which result in collective 5 = 
and S = 1 excitations. Alternatively, the spinons can also 
confine with holons to create charged S' = 1/2 excitations 
with the same quantum numbers as an electron. 

In contrast, in the superconducting state (see Sec- 
tion |^^D] below) the bi are condensed, and then Tif has 
the structure of the Bogoliubov theory of a superconduc- 
tor, with the dis acting as the 5=1/2 Bogoliubov quasi- 
particles; the bi condensate also acts like a Higgs field 
which quenches the A fiuctuations. Indeed, if the bi con- 
densate is strong enough, there is nothing in Eq. (|3.34|) 
which prevents the appearance of gapless, nodal, Bogoli- 
ubov quasiparticles in a d-wave superconducting state 
(see also Section HII D() . In a superconductor which pre- 
serves all lattice symmetries, the nodal quasiparticles will 
appear above a critical doping 5 — 5f as a quartet near 
the center of the Brillouin zone in a strong-pairing to 
weak-pairing transitioniS. On the other hand, the su- 
perconductor could also have broken lattice symmetries 
[i.e. it is a supersolid) and then the transition aX 5 — 5f 
is associated with the pairwise collision of nodal points 
at the reduced Brillouin zone boundariesi^ii^iSS. 

Finally, it is useful to rewrite the expression for the 
physical electron field, Eq. H3.23|l , in the spinon variables: 

b\d- for rji = +1, , . 

t V ( 1 (3.36) 

-bltss'dl^, lor rji = -1 

At low energies, or in the W oo limit, this can be 
rewritten as 



bl+d.s 



f ,d^ 



(3.37) 



since the first and second terms contribute predominantly 
on the A and B sublattices, due to the preferred locations 
of the two hole species. 

How docs this change for the vr flux background? The 
needed form is obtained simply from Eqs. (|3.28() - (|3.34(l 
by shifting according to Eq. H3.14|l . One notes that the 
dis spinons, being minimally coupled to the gauge field, 
indeed propagate in this case according to a 7r-flux hop- 
ping model. Specifically, the shift implies 

/4 



whence 



riiS^e^ 



irjiiT /A 



(3.38) 
(3.39) 

(3.40) 



Uy ^ = r,,<5,e-2'x(«)e+^''--/4^ (3.41) 

for real 5^. Further, it is interesting to note in this case 
that the t' terms in Eq. (|3.30|) vanish with the replace- 
ment Aia = Aia due to destructive interference of the 
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two hopping paths around the plaquette containing the 
two diagonal sites in question. Indeed, one can show from 
symmetry arguments alone that a diagonal hopping term 
is inconsistent with the tt-Hux PSG, using Eq. 13.27|l for 
rotations and Eq. (|3.26l) for inversions. 



C. Connection to slave particle theories 

We comment briefly on the connection of the above dis- 
cussion to the U(l) and SU(2) slave particle approaches 
to the t — J model. The dimer model Hamiltonian above 
can be regarded as a theory of fluctuations around par- 
ticular mean-field states of these approaches. This po- 
tentially allows one to connect to other gauge-theoretic 
scenarios and phases which are not within the purview of 
the dimer model. First consider the U(l) approachiSiiS 
One introduces microscopically holon and spinon fields 
via 

C.S = h\hs, (3.42) 

with the constraint b\h,^ + flj^^ = 1. Eq. should be 

compared with Eq. H3.23|l . Note that here there is only a 
single species of holon on each site, so we must compare 
directly to the dimer model in the large W limit. In the 
path integral Lagrangian, this constraint is implemented 
by a Lagrange multiplier Ai^: 

£o = -E*^'o(6!&. + /I^-l). (3.43) 

i 

The electron hopping {t) and Heisenberg (J) terms of the 
t — J model are decoupled with matrix fields Qij , , Ay- 
leading to quadratic terms in the bosons and fcrmions of 
the form 

A = 5](Q.,5l6^+Q,,ii/,, + A,,iie,,,4,+H.c.) 

(3.44) 

Comparison with the effective Hamiltonians in the 
dimer -|-hole-|-fermion model (for large W) above shows 
that one can obtain the same form at the quadratic 
level by taking a saddle point with iAm — vrji, Qij, Qij 
non-zero only when i,j are on the same sublattice, and 
Aij , Aij non-zero only when i,j are on opposite sub- 
lattices. For the 7r-flux case, Qij, Qij actually are non- 
zero only between sites separated by an even number of 
lattice spacings in the x and y directions (correspond- 
ing to the vanishing of the t' couplings with this back- 
ground, and hence to a four-sublattice structure of the 
bosons). Phase fluctuations of these matrix fields be- 
come the spatial components of the gauge fields of the 
dimer -|-hole-ffermion model. 

Next, consider the SU(2) slave particle 
formulationiiiii^ In this case, one represents the 
electrons by an SU(2) doublet of bosons, bi„, and a 



spin- 1/2 fermion dis, with the electron represented as 

c^s = ^ {bl^d,, + bless'dl) . (3.45) 

Up to a sign change in hi- (due to our trivially different 
conventions), this is the same form as Eq. (|3.37|) . The 
slave particles are constrained by V'LfV'is + blrb- — 0, 
where the doublet 




(3.46) 



and f is the vector of Pauli matrices acting in the doublet 
space. These constraints lead to three Lagrange multi- 
pliers (p= 1,2,3): 

^5^(2) = . Ll^^^^+f,l^A . (3.47) 

ip ^ ^ 

Of particular interest is the p = 3 component of this 
constraint which takes the form 

6l+6,+ - b\_h..^_ + ^ d\A, ^ 1. (3.48) 

s 

Comparing Eq. (|3.48|l with our Eqs. (|3.25|l . (|3.33|l . we see 
that they differ in that Eqs. (|3.25|) . H3.33|l contain an ad- 
ditional fluctuating electric field contribution. This is a 
reflection of the different manner in which the gauge fluc- 
tuations arc realized on the lattice, rather than a funda- 
mental distinction. If we begin with a theory with a con- 
straint as in Eq. (|3.48() , and imagine integrating out some 
high energy degrees of freedom, then the gauge invari- 
ance permits the electric field contribution in Eqs. H3.25|l . 
(|3.33|l to appear in the renormalized theory; thus, the 
main difference is that in our dimer model approach this 
contribution in already included in the bare theory. 

As we will see in Section IIVI in our duality analysis 
the electric field term is a large contribution to the l.h.s. 
of Eq. H3.25|l . which dominates that of the / fcrmions. In 
contrast, in the mean-field saddle points of Refs. Il...l3. 
the dominant contribution to Eq. H3.48|) arises from the 
fermion dj^dj^ term, with the d fcrmions occupying a half- 
filled band (the d fcrmions also contribute a large term 
to the duality analysis of Section IIVI in Eq. H3.33|l , but 
not in a manifestly translationally invariant manner.^^). 
We will comment further on the physical distinction be- 
tween these two approaches in Section IIVI Decoupling 
the hopping and Heisenberg interactions in the Hamilto- 
nian gives quadratic terms for the slave particles of the 
form: 

^smi) ^ Y.c,i;lu,,i;^, + c,blu.,b^, (3.49) 

ij 

with 
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Our dimer+hole+fermion model described by Tigd + 'Hf 
in Eqs. (|2.1|l . (|3.34l) can be described in the above SU(2) 
formulation by taking a saddle point with — — vrji. 
In particular, our physical interpretation requires that we 
take the saddle point value —v ~ W large so that the bier 
holes reside primarily on their respective sublattices. In 
contrast, in previous work^^-J^ on the SU(2) gauge the- 
ory superficially different saddle points were chosen: they 
worked with a spatially uniform value of —iA^Q so that 
only one species of the bosons was preferentially occupied 
in the ground state. However, saddle points in these two 
classes can be mapped onto each other by a SU(2) gauge 
transformation generated by the unitary matrix ir^ act- 
ing only on the rji = —1 sublattice. This gauge transfor- 
mation maps our staggered — i^^o saddle point to a spa- 
tially uniform saddle point. The same gauge transforma- 
tion acting on interchanges the fermion hopping, 
Xij, and pairing, Aij, terms. So to obtain a saddle point 

of consistent with our Eq. H3.34|l we need to take 

the Aij of a form consistent with the Ua in Eq. H3.34|l . 
The Xij "will be proportional to the t times the square of 
the bi condensate. 

When the model is undoped and the holons are gapped, 
{bi) = 0, such states fit into Wen's PSG classification 
scheme for symmetric spin liquid states on the square 
lattice. What are called the zero and tt flux states in this 
paper correspond, in the notation of Ref.0, to the states 
UlCnOOx, and UlCnOnl, respectively (see Eqs (122) and 
(135) of Ref. ,16:). These states are distinct in their 
PSGs from the "staggered flux" and "pi flux" states of 
Refs. IIII12I which are denoted UlCnOln and SU2BriO, 
respectively in Ref. 0. Note that in Wen's discussion, 
the spin liquids UlCnOOa::, and UlCnOnl obey the con- 
straint in Eq. I|3.48|l : consequently the fermion "chemical 
potential" v in Eq. H3.28|) or Eq. H3.34|l is not an arbi- 
trary parameter and must take a value that ensures that 
Eq. (|3.48|) is obeyed. In our discussion, we have the addi- 
tional fluctuating electric field contribution in Eq. I|3.33|) . 
and we treat Eqns. (|3.28l) and H3.34|l as generic effective 
Hamiltonians; large positive values of v, which lead to a 
large spin gap and a small density of / fermions, are also 
allowed. 



D. Superconducting state 

To obtain superconducting states of the doped (or un- 
doped) dimer model, one clearly must condense some 
charged particle. In fact, a conventional superconducting 
state requires a condensation of both flavors of bosons, 
{bier) 7^ for both a ~ ±. Condensation of only one of 
the flavors breaks only the combination of physical U(l) 
charge symmetry (since each holon carries charge -\-e) 
and gauge U(l) symmetry. The orthogonal linear com- 
bination is unbroken in that case. Explicitly, one can see 
this from Eq. H3.37|l and Eq. H3.34|l . When only one bo- 
son is condensed, say (6i+) 7^ 0, then at mean-field level, 
the electron and spinon become equivalent, in this case 



Cis ^ {bi+)dis. Replacing dis by Cis in Eq. H3.34(l and 
dropping the t term which vanishes since one of the two 
bosons is uncondensed, one sees that one has no anoma- 
lous electron pairing terms. One might also imagine con- 
densing not the individual holon fields but instead just 
the quadratic forms (l)\^b-_) ^ and {b\^b'^j_) ^ 0. This 
will certainly break the charge U(l) symmetry, and in- 
deed follows if we condensed both holons individually. 
However, if the individual holons remain uncondensed, 
there is a residual unbroken Z2 gauge symmetry under 
bia Zibia, with Zi = ±1. Thus such a state is an 
anomalous superconductor (SC*), which we do not wish 
to consider here. 

It is simplest to think about the "conventional" super- 
conducting state in the bi and fis variables. Condensa- 
tion of both flavors of holon then simply means that (bi) 
is non-zero on both sublattices. Then, from Eq. H3.23|l . 
we can regard the fis variables as essentially electrons. 
Formally, if the superconducting state does not break 
translational symmetry, the holon condensate will have 
uniform amplitude, 

ih) = \b\e-'^% (3.51) 

so we shift fis — > e~*'^*/is, after which Cis ^ \b\fis- This 
gives 

^? -Hflhs (3.52) 

i 

+ 5]^„e^(^-+^-+^=)/ie.../^,^,,,+H.c., 

where we have dropped the gauge field which is gapped by 
the Higgs mechanism in this phase. The rescaled fermion 
fis can now be regarded as the electron quasiparticle op- 
erator, so that the pairing state can be read directly off 
from Eq. H3.52|l . We see that the pair amplitude obeys 

*L ~'5„e*('^'+'^'+^°). (3.53) 

This is in full agreement with the symmetry discussion 
in Sec lmH 

Note that for the case of a zero flux background, the 
pairing can be of d^2_y2 symmetry. This occurs e.g. for 
a uniform holon condensate (j)i = and Sx — —Sy = 
1 (e2*x(J^) = or for cj), = and 6^ = Sy = 1 

(g2ix(fl) — ly Nevertheless, the spectrum of Tiy is fully 
gapped in the limit relevant to the doped dimer model, 
since there v ^ Sa- This corresponds to a "molecular", 
"BEC" , or "strong pairing" d-wave phase, where it is 
well-known that gapless nodal excitations are absent. 

It is amusing to consider the energetics behind the pair- 
ing symmetry for one example. Let us take for simplicity 
zero background flux, and choose the PSG with Sx = Sy, 
q2ix{R) = 1. To apply Eq. (|3.53l) . we need to understand 
the pattern of holon condensation. For simplicity, let us 
ignore explicit (non-gauge mediated) holon-spinon inter- 
actions, and gauge fluctuations. Then we may consider 
the holons on their own, described by a hopping model. 
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We imagine lowering the chemical potential for holons 
until they condense, the condensate being built out of 
the lowest energy modes. Consider a single holon exci- 
tation. For simplicity, we consider the large W limit, so 
that the flavor 1 holons prefer to sit on the A sublattice 
and the other flavor sits on the B sublattice. For W — oo, 
the lowest one-holon states are degenerate, consisting of 
any superposition of one holon of type 1 on the A sub- 
lattice, or any superposition of one holon of type 2 on 
the B sublattice. One can then apply degenerate pertur- 
bation theory to the holon hopping terms, to arrive at 
an effective model of the type 1 holon hopping on the A 
checkerboard sublattice and the type 2 holon hopping on 
the B checkerboard sublattice. One expect that, allow- 
ing for generalized hoppings in the original model, this 
hopping Hamiltonian takes the form 



Ti-ih — t. 



cff 



(3.54) 



where the prime on the sum indicates that the sum is 
over nearest neighbors of the A(B) sublattice for a ~ 1 
{a = 2). The ellipses indicate other further neighbor 
terms that will generally be present. The single- holon 
eigenstates of Eq. H3.54|l are plane waves, with energy 
Ea{k) = 2tes{cos{kx + ky) + cos{kx — ky)) (momenta 
are defined for the original square lattice). For foff < 0, 



the minimum energy states have either kx 



ki, 



0,7r. 



These two are actually equivalent, and represent states 
with constant bi on the A or B sublattice. For tcff > 0, 
the minimum energy states have k^ = 0, ky = t: or k^ = 
TT, kx = 0, which are also equivalent and represent states 
with bi alternating in sign on nearest-neighbor site of the 
A or B sublattice, e.g. with 6^+ = 1 on sites with Xi,yi 
both even and &i_ = — 1 on sites with Xi , yi both odd. At 
the one particle level, one expects then the condensate 
which forms on doping to be a linear superposition of 
these two flavors, 



(3.55) 



where 6™" is the minimum energy single-holon wave- 
function with support only on the a sublattice, either 
constant or staggered on this sublattice depending upon 
the sign of tes- The nature of the superposition that oc- 
curs - i.e. the coefficients ipi,'ip2 ~ will be determined by 
holon- holon interactions. At a mean- field level, we can 
simply construct a "Landau" effective action for "011 "02 ■ 
Since both flavors of bosons are separately conserved, it 
should have U{1) x U{1) symmetry under rotations of 
each field. Further, it should be symmetric under in- 
terchange of the two flavors, by translation invariance, 
which interchanges the two sublattices. One expects 
therefore an effective (Euclidean) Lagrangian of the form 

-Ccff = -r(|0i|'+|02|')+^^(|0i|' + |V'2|')'-f|0i|'|^2|'. 

(3.56) 

The nature of the condensate depends upon the sign of 
V. For V > 0, one obtains j^il = |'02|, while for v < 0, 



one or the other amplitude vanishes. Since the latter 
case clearly breaks translational symmetry, we focus on 
the former. Then for tes < 0, the phases 4>i are constant 
on all sites of one sublattice, which for tcff > 0, the 0^ 
alternate between and tt on the neighboring (diagonal 
second neighbors on the original square lattice) sites of a 
given sublattice. 

Now let us consider the pairing symmetry. As we as- 
sumed 6y = Sx, we have from Eq. (|3.53() 



(3.57) 



Since i + e^, i + ey are nearest neighbor sites on the same 
sublattice, one finds 



1 for < 0, 
-1 for tcff > 



(3.58) 



Thus the pairing is s-wave or d-wave in these two cases, 
respectively. 



IV. DUALITY ANALYSIS OF THE DOPED 
DIMER MODEL 

This section returns to an analysis the of effective the- 
ory of the total spin 5 = sector presented in Sec- 
tion |n] For this limited theory, we claim that the fol- 
lowing analysis properly accounts for the strong fluc- 
tuations of the compact U(l) gauge field A; a corre- 
sponding, non-perturbative treatment of gauge fluctu- 
ations has not so far been possible in the slave parti- 
cle theoriesiSJiiaJ^JiJ^ The 5=1/2 fermions f^s 
are neglected belowfSi and it is clear that this is surely 
safe as long as the fis remain fully gapped. We noted 
in Section 1111 BI that the structure of our theory allows 
the appearance of gapless nodal fermionic excitations in 
the superconducting phases, and that they do not signifi- 
cantly modify the structure of the resulting superconduc- 
tor. As long as such nodal excitations are present only for 
a, 6 > Sf for which the ground state is a superconductor 
or a supersolid, the following analysis will apply for the 
full range of 6. However, it is also possible that 6f co- 
incides with the position of a superconductor-insulator 
transition: such a critical point will not be described by 
the theories presented below and in 1. 

Before embarking upon our duality analysis, we re- 
mark further on why the neglect of the fis fermions is 
not as dangerous at it might naively seem, and why it is 
amply compensated by the non-perturbative treatment 
of gauge fluctuations. It is useful to refer back to our 
key constraint equations Eq. (|3.25|) . (|3.33|l . and the cor- 
responding constraint equation for the slave particle the- 
ory in Eq. (|3.48|l . In the earlier slave particle theories, 
the dominant, order unity, and translationally invariant 
contribution to the l.h.s. of Eq. H3.48|l is provided by the 
fermion density term d\^d^^. However, by focusing on a 
static gauge field saddle point, these theories do not allow 
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a significant contribution to this constraint equation from 
the l^aEia term on the l.h.s. of Eqs. 1(223, (such 
contributions are only generated upon including fluctua- 
tions about the static gauge saddle points). In our duality 
approach below (and in the previous work^iL^ on the un- 
doped model) a large contribution is given by the A^iJiQ 
term: this is only possible because we account for strong 
gauge fluctuations. The contribution of the / fermionic 
term in Eq. (|3.25() . from either nodal or gapped excita- 
tions, is small (note, however, that the contribution of 
the d fermion term in our Eq. is not smalP^), and 

its neglect appears safe. Physically speaking, the AaEia 
term counts the valence bonds on the links emerging from 
a site, and so our approach here assumes that most of the 
electrons appear in singlet pairs, rather than individually 
in single particle states. 

We proceed with a duality mapping of Hqd as in Ref . 8 
and similar to that described in I. The cosine term in Tit 
is written in the Villain form as 



exp At cos (^Aa(f>icr + crAia 



J2 '^''P 



(4.1) 



'J 1, 



2tI\T 



while that in Ti^ is written as 



exp ( — J At cos ( eapAaAifj 



2JAi 



iBataisAaAiii + inBa . (4.2) 



Here Ba is an integer-valued dual magnetic field on the 
sites of the dual lattice. Then we integrate over the 0io- 
and Aia and obtain the dual partition function 

Z,,= ^ exp(-|^F„2^-z^^F,. 

-^E['^-M-(^^+^^/^w/)v]' 

^ ^ ^ij Jiar J jar J ^ i^^^Jia^ 



(4.3) 



where Jia-f^ — {fia, Jiax, Jiay) represent the integer val- 
ued spacetime currents of the two species of holes, Fa^ — 
{—Ba,—Eiy,Eix,) is the integer valued gauge flux, g — 
1/(JAt) = VAt, and as in I, = tAr = 1/f/Ar. The 
iirFar term will be innocuous in the present subsection, 
and we will return to a consideration of its effects in Sec- 
tion llV Cl when we examine the signature of the electronic 
pairing. 



• i, • * I * * 



1/4 1/4 
••••< 

X k'^ X 3/4 1/2 3/4 1/2 



X X 



(a) X,, 



1/4 1/4 
(b)X 



+1/8 -1/8 +1/8 -1/8 
XXX 

-1/8 +1/8 -1/8 +1/8 
XXX 

+1/8 -1/8 +1/8 -1/8 



(c)Z 



■in 



FIG. 3: The values of the only non-zero components of the 
fixed field A'a/j, Va, and Zi^. The circles (crosses) are the sites 
of the direct (dual) lattice. In (c), only the fi = t component 



We solve the constraints in Eq. (|4.3|) by parameterizing 

Jiafi — ^fiuxAi/ho^fjX 

Fat, = At,Na+ Xat^-^Crbaap, (4.4) 

a 

The integer- valued fixed field Xaf, is defined as in Ref. |^ 
and is shown in Fig 13 it obeys the constraint 



(4.5) 



The degrees of freedom are integers 600-;^ on the links of 
the dual lattice and integers Na on the sites of the dual 
lattice. 

Now we promote these integer-valued fields to real 
fields by introducing a vortex fugacity yy and a monopole 
fugacity i/m- Then 

ac |_ 

- ZTT ^ \^rNa + Xar <jh 

ia 

+ Vv ^ C0s(A^l9a<T - 27r6aCT^) 
acr 

a \ a / . 

This is a theory of 2 vortex fields, 7/1+ = e'''+ and 
-!/)_ = e^^- which are coupled to 2 non-compact U(l) 
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gauge fields 6+^ and There is also a Higgs field Na 
which breaks the staggered component of the U(l) gauge 
symmetries, as will become more explicit below. 

We will henceforth not carry through the off-site Ay- 
terms explicitly in our analysis. As in I, we will assume 
that their effects can be absorbed into a renormalization 
of the value of e^. This is certainly the case for short- 
range Ay , while for long-range Coulomb interactions, 
will acquire momentum dependence which was noted in 
I. Our numerical results below will be restricted to the 
short-range case, but we do not expect significant modi- 
fications for the long-range case. 

To allow us to make direct contact with previous work 
on the undoped dimcr model, we use the parameteriza- 
tion of Xafi in Ref . |8i in terms of curl- free and divergence- 
free fields (shown in FiglSJ 



in its final form 



iX- 



(4.7) 



Inserting this in Eq. H4.6|l and shifting Na Na — 3^a we 
obtain 

Zqd = n y" d-^aati J d^aa J dNa CXp 
i(7 

+ yv^ cos(A^i9 ^Trbaafj.) 



J2 COS ( 27rNa - J2 '^^'^'^ - ^^y^ 



(4.8) 



where £a^ = e^^\A,yZtx is given by 



f 



(0,1,-1). 



(4.9) 



It is convenient to introduce uniform (A^) and 'staggered' 
{Cfj,) gauge fields defined by 



^a — fi ~ Aafi Cafj.. 



(4.10) 



The term proportional to g in Eq. 1)4. 8|l effectively breaks 
the 'staggered' gauge symmetry associated with under 
which ip± have opposite gauge charges. The field Na is 
the phase of the Higgs field which breaks this symmetry, 
and as is conventional, we use a gauge transformation to 
set this field equal to 0. Also, we transform from the 
hard-spin to soft-spin variables to obtain the dual action 



a 

^ ^ [e^iyxAiyCaX - r}iW5^rf 

i 

^ f X] ^ ^afi)'^ + 2z7r^ Car 

a a 

aa a 

+ yE[v^:+.,<^e^^^'^°^+-^-)^..-t-c.c 



(4.11) 



The similarity between Zgd above, and the dual theory 
of bosons on a square lattice in I should now be evident. 
The latter theory had a single vortex species, tpa, coupled 
to a single non-compact U(l) gauge field Aa^. Here in 
Eq. (|4.11|) we have two vortex species, 'ipa+ and V'o-, 
but they are coupled together by the monopole fugacity 
term proportional to y™, effectively reducing the theory 
to that of a single vortex field. The present partition 
function also has 2 non-compact U(l) gauge fields, 
and Ca^, but the Ca^ field has acquired a mass from the 
Higgs phenomenon, as is clear from the term proportional 
to g. This massive gauge field also has a complex Berry 
phase term (with 2i7r in Eq. (|4.11|l '). which is unimportant 
here because the Cafj. fiuctuations are quenched. The 
physical meaning of this quenching was already discussed 
at the end of Sec. ^1 i^i this representation, the regime 
of interest where the dimer resonance move has relatively 
small amplitude corresponds to large g, and the tt gauge 
fiux favored by this term (responsible for the complex 
Berry phase term just mentioned) is expelled. 

Given the close similarity between Zd in I and Zqd in 
Eq. H4.11|l . it is evident that low energy fluctuations in 
the vicinity of transitions between superfluid, supersolid, 
and insulating phases of the quantum dimer model are 
described by the same continuum quantum field theo- 
ries as those in I. The most important fact determining 
the character of these field theories is the particular pro- 
jective representation of the square lattice space group 
(PSG) that is realized by the saddle point under con- 
sideration. In particular, all we need is the value of the 
effective dual 'magnetic' fiux /, and the corresponding 
value of the unimodular complex number uj in I. Deter- 
mining this requires a careful symmetry analysis of the 
appropriate saddle point of Eq. H4.11|l . and this is carried 
out in the following subsection. 

Before turning to this symmetry analysis, it is useful 
to make explicit contact between the dual partition func- 
tion in Eq. (|4.11|) . and previous dual representations^'® of 
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the undoped dimer model. The hole density vanishes in 
this limit, and so the flux is zero; this is achieved by 
setting = in Eq. (|4.1HI . Furthermore, the insulating 
behavior requires a strong condensate of the il'a± flelds, 
and we can focus on the phase fluctuations of the con- 
densate by writing i/)a± = e'^°± . We also set the massive 
fleld Cap — 0. Then the action Zqd maps at low energies 
to 



qd 



n / ^"^"/^ n / '^^'^'^ '^^p 

ail aa 



yv_ 

2 

Vv 



afj, 

cos (x 



(4.12) 



The flrst two terms show that the gauge field is Hig- 
gsed by Xa+ + Xa-, and so is innocuous. The last two 
terms constitute a sine-Gordon model for Xa+ — Xa-, with 
on-site offset ya'. this is precisely the dual representation 
of the undoped dimer model found earlieri2i£. 



A. Symmetries 

As in I, we will analyze the symmetry properties of the 
dual dimer model theory in Eq. (|4.11() at a commensurate 
density of holes appropriate to a proximate Mott insula- 
tor. We assume this has the rational value 5/ which we 
parameterize as 



5i 
2 



P 

7 

9 



(4.13) 



where p and q are relatively prime integers. In general (as 
in I), the density of holes, J, in the theory in Eq. (|4.1HI 
is determined by the parameter 7? , and in the superfluid 
or supersolid phases we may have S Si. 

Proceeding as in I, we set the gauge field equal 
to a saddle point value Aafj. such that the flux is equal 
to the density of each species of hole (+ or — ). So we 
choose Aar — Aax = and 



A -^a 



(4.14) 



For the staggered gauge fleld we have the saddle point 
value Cati = l^aM with A = 1/(1 -t- A/e^g). 

The main new subtlety here (beyond that of I) in the 
symmetry analysis of Eq. H4.11|l is the presence of the 
flxed background flelds Eafj. and J^q. These fields are not 
explicitly invariant under the square lattice space group, 



and so require additional transformations of the vortex 
flelds upon operations of the space group. These ad- 
ditional transformation will modify the needed PSG, as 
shown below. 

Carrying out the space group operations in the pres- 
ence of the Aa^, San, and it can be shown that the 
action remains invariant under the following fleld trans- 
formations, which are the analog of the relations in Sec- 
tion II.B of I: 



Ty : ^_(fc,,fc^)-> V+(fc,,fc^)e-^'=''e-'-/4 
Ta; : ip+{kx,ky) ^ ij-{koi:,ky - TrSi)e 



t/4 



Tx : ^P.ikx,ky)~.^+ikx,ky~ 7:61)6-''''' e'^/* 
R-k/2 ■ i^+{kx,ky) 

1 

- ^^(ky + TTnSi,-kx + TTmSi)Lu ™" 
q ^-^ 

m,n— 

i?^/2 : ij-{kx,ky) (4.15) 
1 

— ■^lJ-^-{ky + Trn5j,—kx + Trm6i)LLj ™" 

^ m.n—Q 

To understand the degeneracy of the spectrum associated 
with these transformations, it is useful to transform to a 
basis of states which are eigenvectors of Ty. These are 
expressed by the flelds 



V2 
v2 ^ 



(4.16) 



Then 



y . Tpi{kx,ky) ^ iji{kx,ky)e 



Ty „, 

Ty : iJ2{kx,ky) ^ ~ip2{kx,ky)e-'''y 
Tx : Mkx,ky)^Mkx,ky-7Tdi)e-'''''e'-^^ 
Tx : Mkx,ky)^Mkx,ky-7TSi)e-'''''e-'^^^ 
1 

R : Mkx,ky)^- y to-"'" 

m,n— 

+ i^2\{ky + imSj, —kx + TTmdi) 

1 

R : Mkx,ky)^— V (4.17) 
q ^-^ 

m.n—O 

[tp2 + i'tpi] {ky + nnSi ,—kx + ttuiSi ) 

Now we observe that the Ty eigenvalues of the tpi flelds 
at momentum k are equal to those of the '02 flelds at 
k + (0, tt). So these flelds will mix with each other, and 
cannot lead to orthogonal eigenmodes. It is convenient to 
introduce a new collective fleld \E'(fc) which equals tpi at k 
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and ■02 at + (0, tt), and which is identified by the eigen- 
value e"*'"!' of Ty. Then, ignoring unimportant phase 
factors, we have under the action of T^'. 



(4.18) 



From this equation we see that the PSG of the present 
fields has the same structure as that in Section I of I 



with the key parameter / specified by 



/ 



P 

q 



(1 



p 

7 

q 



(4.19) 



(4.20) 



where p and q are also relatively prime integers. Note 
that the value of / is half the density of electrons i.e. it 
is the density of Cooper pairs in the Mott insulator. So 
the low energy field theory of the present dimer model 
is identical to that presented in I but with the PSG de- 
termined by the density of Cooper pairs. We have also 
explicitly diagonalized the spectrum of ipaa vortices in 
Zqd in Eq. (|4.11|) and verified that it did indeed have the 
degeneracy required by such a PSG. 



B. Mean field theory 

We have seen above that the low energy continuum 
theory of the quantum dimer model is identical to that al- 
ready presented in I. We have already presented a mean- 
field analysis of this continuum theory in Section II. D of 
I, and so do not need to repeat it here. Rather, we will 
work directly with the lattice action in Eq. (|4.11|) and 
follow its mean-field phase diagram with increasing hole 
density, S. 

To do a mean-field analysis of our dual dimer model 
we consider the free energy corresponding to Eq. H4.11|l . 
Rescaling the fields ipaa, and neglecting an overall con- 
stant, we can choose u — yy — I and obtain the free 
energy (with r = s — 2) 



-2i:i{Aa^+(jCa 



E 



1, 



aa '- 

a 

a 
a 



vY^\^a+?\i>a-\ 



(4.21) 



For simplicity we will also set ym equal to 1, and, as 
discussed below Eq. p. 8(1 . we can safely choose W = Q. 
We will set v equal to zero and do not discuss states 
which break the a ^ —a symmetry. Such states break 
the vortex-anti- vortex symmetry on average, and emerge 
only for sufficiently large v. 

Wc note that we have dropped the 2iTT Car term 
in Eq. (|4.11|) , which was a consequence of the sign J > 
in the dimer 'resonance' term Eq. (|2.9() . The massive 
Car field has little effect on our analysis below, and it 
therefore appears reasonable to ignore this term in our 
me an- field theory. 

This leaves us with the hole chemical potential (the 
'magnetic field') H (in the superfiuid phase, this picks a 
hole density 5 — 2H), the coupling constants and g, 
and the control parameter r. For the phase transition 
between the superfiuid and the insulating phase to be of 
second order we have to choose small. In our units, 
we take = 0.04. 

We can now map out the phase diagram for zero dop- 
ing. Instead of using r and g as parameters, let us use r 
and the already introduced (below Eq. H4.14|l ) parameter 



A = 



1 



1 + 4/e2g 



(4.22) 



Note that the range < g < oo gets mapped onto < 
A < 1. Also, in the superfiuid phase, and close to the 
second order phase transition to this phase, the saddle 
point approximation for the massive field Ca^ becomes 

exact such that iJstaggered = erMi/A^Cay = ±A/2. 

As can be seen in Fig01 for large enough r we are al- 
ways in the superfiuid phase, which is characterized by 



the vortex vacuum {ipafj 



0. As we decrease r, we 



enter an insulating state with {4'an) 7^ via a second or- 
der phase transition. Depending on the value of A, this 
can either be a columnar valence bond solid (VBS) state 
(B), or a state with site density wave order (CDW) at 
wavevector (7r,7r) (^4). It turns out that the plaquette 
state (C) always has a free energy slightly larger than 
that of the columnar VBS state, for the parameters cho- 
sen here. All these results are consistent with the mean 
field phase diagram of our field theory for g = 2 in I. 

In Fig 21 (and Figs I6I9|I we have indicated the links of 
the direct lattice by solid lines. The size of the (black) 
bars living on these links is proportional to the vortex 
kinetic energy 



(4.23) 



which we can physically interpret as the mean spin ex- 
change energy of the electrons on the direct lattice. The 
mean electron ring exchange energy determined by 



(4.24) 



is depicted by (blue) boxes whose size is chosen to be 
proportional to this quantity. If the exchange energy is 
negative we have used a crossed box. Finally the (mean) 



0.3 



0.1 



-0.1 



B 



favored 
over 




FIG. 4: Mean field phase diagram for the doped dimer model 
with zero doping. For large r our theory favors a vortex vac- 
uum such that we are in the superfiuid phase. Decreasing r 
we enter one of the (condensed) insulating states. While for 
smaller A we find a columnar bond ordered state [B) whose 
free energy is slightly lower than that of a plaquette state (C), 
for larger values of A we get a CDW state {A). The lines in- 
dicate the links of the direct lattice. The black bars on these 
links are proportional to the spin exchange energy or the elec- 
tron pairing amplitude (see text). Similarly, the size of the 
(blue) boxes depicts the ring exchange amplitudes which are 
located on the dual lattice. If this amplitude is negative we 
have used a (tiny) crossed box. Finally, the state A also shows 
a finite mean hole density which is positive on one sublattice 
(red bullets) and negative on the other (crossed circles). As 
discussed in the text this CDW state is unphysical, and we 
exclude the corresponding parameter values. 



density of holes on the sites of the direct lattice is pro- 
portional to 



A A 



(4.25) 



For H = (zero doping) we expect this quantity to van- 
ish, which indeed is the case for the bond ordered or 
the plaquette state. However, the CDW state in Fig^ 
shows a checkerboard pattern for Bi, which takes oppo- 
site values on the two sublattices, and are depicted by 
(red) bullets or crossed circles depending on whether Bi 
is positive or negative. Clearly this CDW state is not in a 
physically interesting regime of couplings for the original 
dimer model as it applies to the cuprates (the electron 
density is uniform in the Mott insulator at half-filling) ; 
the lesson to be learned from this is that in our model 
we cannot choose g too large. For concreteness, let us 
choose g = 2/e^ = 50 such that for sufficiently small r 
we are always in the columnar VBS state. 

To discuss doping holes into an insulator, we will start 
from a columnar bond ordered state and follow its fate as 
we increase the chemical potential H. First of all, it is in- 
structive to determine the critical value of r at which we 
have the phase transition from a superfiuid to an insulat- 
ing state as a function of H. This transition occurs when 
the tpan condense, and so a lower bound on the critical 




FIG. 5: Lower bound for the superfiuid to insulator transi- 
tion, Vc as a function of the hole density, S, in the superfiuid 
state {5 — 2H). There is a second order transition at the 
points shown into an insulator with the same density. For ir- 
rational 5, this insulator is a fioating Wigner solid. The ratio- 
nal 5 Mott insulators are expected to survive in fixed density 
"Mott lobes" extending over a finite range of the "chemical 
potential" H, and these can pre-empt the transition into the 
incommensurate fioating solids. The re-entrance of the curve 
above at large // is in a regime of such large 5 that we do not 
expect our small doping theory to remain valid. 



value of r = Tc is obtained by locating the point at which 
the lowest eigenvalue of the quadratic part of the ijja^ 
free energy in Eq. (|4.21(l turns negative. A graph of the 
lower bound of rc{H) is shown in FigEl Determination of 
this lowest eigenvalue requires the Hofstadter spectrum 
of i/ja/i , and so the graph in FigO has a shape similar to 
the ground state energy of the Hofstadter butterfly. For 
incommensurate values of H , the transition at the points 
in Fig [SI is into a floating Wigner solid with hole density 
5 = 2H. We expect that this will often be pre-empted 
by a flrst order transition into a nearby commensurate 
Mott insulator with hole density S ^ 2H , but instead 
at the rational value determined by the Mott insulator. 
In the dual language of the 'superconductor', this Mott 
insulator has expelled some of the 'magnetic flux' H . 

In the following we will choose r — —0.1 such that we 
are in the insulating state at zero doping, and enter the 
superfluid phase at around 5 « 0.15. Our mean field 
theory also predicts that at high enough doping concen- 
tration we are back in an insulating state, but we do not 
expect our duality analysis to be valid at such large 5. 

We do not make any pretense at completeness in the 
following phase diagrams: they are merely representative 
structures obtained for a sample set of parameters, and 
other parameters can clearly give a wide variety of other 
phases. However, all such phases will obey the general 
constraints we have delineated in I. 

In Figs. I6I9I we show the states for doping 5 = 
0, 1/32, 1/16, and 1/8. While at zero doping we are in a 
columnar VBS state with vanishing hole density on the 
sites, at doping 5 = 1/32 we find a modified plaquette 



17 




FIG. 6: Insulating state with columnar valence bond solid 
(VBS) order at zero doping. This state is obtained by a 
mean-field minimization of the free energy in Eq. (14.2111 with 
parameter values as discussed in the text. The significance of 
the symbols on the square lattice above (and in Figs 1719^ are 
as in Fig 2] 

state and which shows a concentration of the mean hole 
density at the corners of one plaquette of its 8 x 8 unit 
cell. At the center of this plaquette frustration leads 
to (V'a+) = (V'a-) — 0. The variation of charge den- 
sity in the unit cell clearly results in a 2-dimensional 
charge density wave with wave- vectors k^^ = (27r/8a,0) 
and kj, = (0, 27r/8a) where a is the lattice spacing. It 
should also be noted that as for the plaquette state (C) 
in Fig21 the state depicted in Fig0is invariant under ro- 
tations by 90 degrees about the plaquette with (V'acr) = 0. 

Doubling the doping concentration S, we find the state 
depicted in Fig|Sl which is similar to the previous state 
but has a unit cell of half its size and leads to a mod- 
ulation of the charge density with wave-vectors k± = 
(27r/8a, ±27r/8a). Finally at doping S — 1/8, we find the 
state depicted in Fig |51 which is also invariant under ro- 
tations by 90 degrees about a plaquette with (ipaa) — 
and has a 4 x 4 unit cell. Since we are now close to the 
second order phase transition to the superfluid phase, 
the mean hole density is now distributed essentially uni- 
formly over the lattice. However, there is still a weak 2- 
dimensional hole density modulation with wave-vectors 
k^ = (27r/4a,0) and ky = (0,27r/4a) 



C. Pairing symmetry 

The analysis of this section has been carried out al- 
most entirely in a dual representation of the underlying 
electronic degrees of freedom. While this has the benefit 




FIG. 7: An insulating state with a 8 x 8 unit cell for doping 
S — 1/32. The parameters are the same as for the 5 = 
state in Fig |S| apart from the change in H, the hole chem- 
ical potential. This state is invariant under rotations by 90 
degrees around the plaquette with {"il^acr) = 0. The mean hole 
density indicated by the red bullets is concentrated near this 
plaquette. 

of properly accounting for the non-perturbative dynam- 
ics of the compact U(l) gauge theory Hqd in Eq. (|2.1|l . 
it does obscure the physical interpretation of the phases 
of the theory. We have already shown in I how the den- 
sity wave order can be thoroughly characterized by the 
dual vortex degrees of freedom. In this subsection we 
want to determine and analyze the corresponding dual 
representation of correlators which characterize the su- 
perconducting order. 

We already discussed the issue of the pairing symme- 
try in Section IIII Al and found that certain key deter- 
mining factors could not be specified within the physics 
of the dimer model alone. In Eq. H3.7|l we related the 
electron pairing field to the product of microscopically 
determined co-efficients fia and degrees of freedom of the 
quantum dimer model. Let us consider the two-point 
correlation function for the pair field, at two pairs of 
nearest-neighbor points 12, 34, where 2 is displaced in 
direction a from 1, and 4 is displaced in direction (3 from 
3 (see Fig EH for an example). Then we can write the 
pairing correlator as 

(^12*34) = Ma/^/3Gl2;34 (4.26) 

where Gi2;34 is a two-point correlator of the quantum 
dimer model given by 

Gi2;34 - {ClCu) , (4.27) 
with the gauge-invariant hole pair annihilation operator 

C12 = e'^^^^^^e-'^'+e-'^^-. (4.28) 
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FIG. 8: Insulating state evolving out of state in Fig |7| for 
doping S = 1/16 (other parameters as in Fig El- As for the 
state depicted in Fig[7|this state is symmetric under rotations 
by 90 degrees. However, the unit cell is now 4x8. 




FIG. 9: Insulating state for doping 5 = 1/8 with a 4 x 4 unit 
cell, with parameters as in Figs. I6I8I Now we are close to the 
phase transition to the superfluid phase and the mean hole 
density is more evenly distributed. 



We choose sites 1,3 on the A sublattice (77^ — +1) and 
sites 2,4 on the B sublattice {rji = —1) for concreteness. 

It now remains to evaluate G'i2;34 within the context of 
duaUty analysis of Section llVl We will show that this cor- 
relator factorizes in the limit of large separation between 
the points 12 and 34, with no dependence on the rela- 



FIG. 10: Computation of the electron pair correlation func- 
tion. The sites 1,3 have rji — +1, while 2,4 have Tji — — 1. 
The current Jai^t is equal to unity along the line connecting 3 
to 1, and is zero elsewhere; other currents between two points 
are defined similarly. We also define loop currents e.g. Ji243,j 
is unity only along the loop around the points 1,2,4,3. An 
electron singlet pair is created on the link 1,2 and annihilated 
on the link 3,4. 



tive orientations or centers of 12 and 34. This is indica- 
tive of superconducting Off-Diagonal Long Range Order 
(ODLRO). The symmetry of the pairing state is then de- 
termined by the dependence of (^'i2^'34) on either pair 
of neighboring points 34 or 12, separately. Consequently 
the pairing symmetry of the superfluid ground state is 
determined by the /i^ factors controlling the particular 
PSG realized by the dimer model, and by the depen- 
dence of G'i2;34 on 12 or 34 separately. We will see that, 
for the large g limit on which we focus, G'i2;34 is also in- 
dependence of the orientation or position of either pair 
of points when they are well-separated. Thus the pairing 
symmetry is determined entirely by the fia factors. 

The correlator Gi2;34 can be computed in our dual rep- 
resentation by inserting the operator in Eq. (|4.28l) . and 
the corresponding operator at sites 3,4, into the trans- 
formations leading to Eq. I|4.3|l . After integrating over 
the degrees of freedom on the direct lattice as before, the 
correlation function is expressed as 

Gl2;34 = ^ exp I - 1 ^ F^^ -iT^Y^ Far 

i 

/ icr 

afi \ (T / 

(4.29) 

Here the currents Ji2^ etc. are fixed background cur- 
rents defined in Fig ^] associated with the presence of 
the source terms of the correlator. Notice that the delta 
function constraints in Eq. H4.29|l differ importantly from 
those in Eq. H4.3|l ; consequently the solution of these con- 
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straints in Eq. H4.4|l is replaced by 

Fa^ = A^Na+Ba^ + Xa^-Y^abaa^., (4.30) 

a 

where Bap. is a fixed background integer-valued field 
which obeys 



£fj.u\A,yBa\ — — Jl243/J 



(4.31) 



Therefore, is the "magnetic" field associated with the 
current loop in FiglTUb: a convenient choice is to take only 
Bar non-zero along the links which pierce the current loop 
in Fig llOb. Including the additional offset Sa^, we can 
continue the analysis of Eq. H4.29|l as before, and obtain 
an expression corresponding to that in Eq. (|4.6|l . 

We now split the baafi gauge fields into its uniform and 
staggered components as in Eq. (|4.1U|) . Manipulating the 
expressions in the action in Eq. H4.29|l wc find that they 
can be split into those involving the uniform (Aa^) and 
staggered (Ca^) components respectively. While the re- 
sulting action has decoupled contributions from Aa^ and 
Cap, it should be kept in mind that their fluctuations are 
not truly independent. In particular, because the batrp 
must be integers, the Aap and Cap must both be either 
integers or half- integers. As we argue below, the Aap 
field fluctuates strongly in the superfluid phase (so that 
it can effectively be considered a continuous real field), 
and so it appears reasonable to ignore such a constraint. 
However it is possible we are overlooking some subtlety 
here, and it may be worthwhile to re-examine this issue 
in future work. 

For the uniform component wc obtain the terms in the 
action 



Sa[J] ~ ^ {2e^,^\A^Aa\ + Jsip + Ja2p 



2HS„ 



(4.32) 

Notice that the source term here is the sum of the two 
currents fiowing from one pair field to the other, and so 
upon coarse-graining it looks like the current of a charge 
2 object {i.e. a Cooper pair) moving between 12 and 
34. In the superfiuid phase, the uniform gauge field Aap 
is in its Coulomb phase because the vortex fields ijjaa 
are gapped. So we can easily evaluate the expectation 
value of the action in Eq. I|4.32|l by treating Aap as a real 
variable controlled by the usual Maxwell action. In this 
manner we deduce that the contribution of the uniform 
component of the gauge field to G1234 is 



JVAape-^-^J^ 



fDAap 



-5a[0] 



(4.33) 



exp ( (A^(J3i^ + J42p)) _^^2/^2 i^^J^l'^ + -^421.)) 



The divergence of the currents is only non-zero at the lo- 
cation of the Cooper pair source terms, and so we see that 
the correlator is the magnetic energy of two 'monopole' 



sources at these locations. The square root of this corrre- 
lator in the limit of infinite separation is the superfiuid 
order parameter, and this is given by the finite self-energy 
of a single monopole. In addition, there is a 1/r inter- 
action between the monopoles, and this power-law decay 
arises from the superfiow fluctuations in the dual super- 
fluid. 

We turn next to the contribution of the staggered field 
Cap to the action in Eq. H4.29|l . In the gauge Na — 
the relevant contribution to the action is 

Sc[J] = \Y. (^-M + ^ap - 2Capf 

a 

+ ^ (Bar + Xar - 2C ar) (4.34) 

a 

+ X! {^^l^vxAyCaX + J-i\p + J2V - ^.-q^W^^rf ■ 
a 

We now shift 2Cap by Bap + Xap, and using Eqs. H4.5() 
and H4.31|l the action above can be written as 

Sc[J] = 2gJ2cl~2znY,Car 

a a 
a 

-r,,{2W + l)Spry . (4.35) 

The most important property of this expression is that it 
only involves the local currents ■Ji2p and -Jisp which are 
entirely localized at the Cooper pair source terms. So 
when these sources move far apart, the correlator factor- 
izes into two contributions, one for each order parame- 
ter insertion. Furthermore, because the staggered gauge 
correlations of Eq. (|4.11|) decay exponentially with sep- 
aration, it follows that these insertions can be evaluated 
independently of each other. As discussed earlier, this is 
indicative of ODLRO. 

Next consider the dependence of G'i2;34 upon 34, keep- 
ing 12 fixed and far away. For two sets of points 34 and 
3'4', one has 

(4.36) 



G 



12;34 



G 



12;3'4' 



-Sc[Jl2,J; 



3'4'l 



For large g, one clearly sees that the numerator and de- 
nominator are dominated by Cap = 0, independently of 
J34 and J3'4'. Hence there is no directional or other de- 
pendence of G'i2;34 on either pair of points, and we can 
conclude that Gi2;34 reaches a featureless constant in the 
superfiuid phase as 12 and 34 move apart from each other. 
The pairing signature in Eq. H4.26|l is determined by the 
/ifj, as promised. 



V. COLLECTIVE S = 1 'TRIPLON' 
EXCITATIONS 

We discussed the inclusion of non-zero spin excitations 
in the dimer model in Section IIIII but focused entirely 
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on fermionic S* = 1/2 degrees of freedom. We also noted 
in Section 1111 B\ that neutral 5* = 1/2 excitations did not 
survive in the insulator: instead they confined in pairs to 
create stable, bosonic 5 = 1 'triplon' excitations. It is to 
this triplon excitation that we now turn our attention to. 
In principle, the triplon can also survive as a stable exci- 
tation into the supersolid or superfluid phases, provided 
energy and momentum conservation constraints prevent 
its decay into pairs of S* = 1/2 Bogoliubov quasiparticles. 

Microscopically, we can also see the necessity of includ- 
ing the triplon excitation by referring to previous anal- 
yses using bond operatorsi^iSSiS^. As discussed therein, 
two electrons on sites at the ends of a square lattice link, 
in addition to forming a spin singlet valence bond (a 
dimer), can also be in a higher energy S — 1 combi- 
nation. It is the motion of this S ~ I state that we wish 
to examine here. 

Clearly, establishing the stability of the triplon requires 
an accounting of fluctuations of the compact U(l) gauge 
field A. So we need to work with dual fields of Section lTVI 
to account for the 5* = sector. It is cumbersome to 
carry triplon through this duality analysis, and so we 
will attempt to guess its effective action using symmetry 
considerations. 

Our approach will be to combine the approach of 
Ref. with the results of I and of Section IIVI As in 
Ref . l24i we represent the triplon by a real vector field fii , 
where the arrow denotes a vector in spin space. The spin 
operator on site j is related to nj by 



fields we can also construct the density operators 



5j cx rjiUj 



(5.1) 



so, on its own, fii is a measure of spin correlations at 
the commensurate antiferromagnetic wavevector (7r,7r). 
Given the short range antiferromagnetic couplings be- 
tween the spins (or, alternatively, using the results the 
microscopic calculation of Ref. , we follow Ref. |24' in 
postulating the following phenomenological quantum lat- 
tice model for the Hi: 



9-1 



f=0 



(5.3) 



where to, n are integers, representing Fourier components 
of 'density'-like observables (such as pairing amplitude, 
exchange energy, or site charge density) at wavevector 
27r/(TO,n). Finally, from the we can Fourier trans- 
form to real space and obtain 8p{v\ the value the 'den- 
sity' at the sites, links, and plaquettes of the square lat- 
tice; the manner in which this is done, and the conven- 
tions for r, are described below and in Eq. (2.26) of I. 

We can now compute the effect of the 5 = sector 
on the triplon n excitation just as in Ref. |3- We as- 
sume that the variations in bp{r) modulate the exchange 
constants between (and the amplitudes of) the n^, and 
so write down the simplest local couplings between these 
degrees of freedom: 



Sn 



dr 



i 



\2 Sp{r 



+ A 



3 2^ '5p(rj 



(5.4) 



We can now analyze Sno as in Refs.Hmil For 
static condensed ipi^ as obtains in a commensurate Mott 
insulator, we have a corresponding static modulation in 
5p{r). Its influence on the triplon excitation spectrum 
can be as computed in Ref. |2j[ As was shown there, sim- 
ple and natural choices for the modulations agree with 
neutron scattering observations2SiSS*2L2£ over a wide en- 
ergy range. We leave for future work the extension of 
these results to the case where the (p^ are fluctuating. 



VI. CONCLUSIONS 



5. 



nO 



in') 



dr 



u 
4 



(5.2) 



Here s is a parameter which determines the gap towards 
triplon excitations, u is a quartic non-linearity, and c is 
a velocity. In the absence of any coupling to the 5 = 
charged excitations, Eqs. H5.1|l . H5.2|l predict that the the 
lowest energy triplon excitation is at wavevector (tt, tt). 

It is now useful to characterize universal aspects of the 
vortex theory of Section lTVI using the methods of I. From I 
we learn that ioi f = p/q {p, q relatively prime integers), 
the low energy physics of the 5 = sector is captured by 
q vortex fields (/?£,£ = 0, 1 ... (q — 1). The effective action 
for the (pf fields is as described in I. From these vortex 



The primary purpose of this paper was to present a 
complete treatment of the interplay between VBS order 
and superconductivity in a model of a doped 5=1/2 
quantum antiferromagnet on the square lattice. We 
aimed to do this in an approach which properly accounted 
for strong gauge fluctuations even at non-zero hole con- 
centrations, S. While numerous previous studies have 
studied the properties of effective gauge theories of quan- 
tum antiferromagnets, essentially all have treated gauge 
fluctuations only at (5 = (see however Ref. 29). As 
we recall briefly below, several key characteristics of our 
theory were special to 5 = 0. 

Our analysis here was carried out in the context of the 
doped quantum dimer modeli, which served as a conve- 
nient effective theory for the spin 5 = sector. By a 
duality analysis of this model in Section HVl we obtained 
a dual theory of vortices in the local paired-electron su- 
perconducting order. As in a previous paperi (referred to 
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as I) on pure boson models, we found that these vortices 
appeared in g > 1 'flavors' determined by the electron- 
pair density in a proximate Mott insulator with 6^0. 
There is much evidence in the cuprates that they are 
proximate to a Mott insulator with hole density 6 — 1/8, 
and by Eq. (|4.2()(l this translates into a value q — 16. 
With the knowledge of q, we can then freely borrow over 
essentially all of the results of I on the fluctuations of su- 
perflow and charge/VBS order across superconducting, 
supersolid, and insulating phases. The q vortices trans- 
form under projective representation of the square lat- 
tice space group (PSG) which strongly constrains their 
effective action. Any impurity-induced pinning of a vor- 
tex also induces VBS/charge order by breaking the PSG 
among the q vortices. We presented a number of results 
on features of impurity induced pinning of vortices in I, 
and the present paper has provided the promised justi- 
fication for their application to scanning tunnelling mi- 
croscopy (STM) experiments on the cuprates. "'"■''^•'^^■'''^■■^^ 

We also showed how the framework of the quantum 
dimer model could be extended to include fermionic and 
non-zero spin excitations. This led us to an unexpected, 
and fairly explicit, connection between the dimer model 
and previousM2Jii2, slave-particle U(l) and SU(2) gauge 
theories of quantum antiferromagnets. This connection 
was established in the context of mean-field saddle points 
in which all gauge fluctuations were quenched. Our 
present analysis makes it quite clear that neglect of such 
gauge fluctuations is dangerous at least at small 6: all 
the PSG properties of vortices described in the previ- 
ous paragraph emerged upon a careful treatment of such 
fluctuations at J ^ 0. 

The relationship between these different approaches to 
quantum antiferromagnets is neatly highlighted by the 
key constraint equation in Eq. H3.33|l 

A„^.„ -I- blb^^ - blb^_ + ^ dld^^ = 1, (6.1) 

s 

obeyed at every site of the lattice. Here Eia is the integer- 
valued electric field of a compact U(l) gauge field which 
is a measure of the number of singlet valence bonds on 
a link of the direct square lattice. The bi± are bosons 
representing holes which carry opposite gauge charges, 
and the d]^ are S = 1/2 fermionic spinon degrees of 
freedom. (It must be emphasized that this is merely a 
kinematic description of theory, and the ultimate quan- 
tum numbers of the emergent quasiparticles may well 
be diffe rent.) In the SU(2) gauge-theoretic approaches 
of Refs the first valence bond term in Eq. H6.1() 

is not explicitly included, although it will be generated 
upon renormalizing the theory to a lower energy scale; 
the constraint is satisfied primarily by a half-filled band 
of d fermions. In our present dimer model approach, in 



contrast, the fluctuating electric fleld term plays a cen- 
tral role in Eq. 1)6.1(1 : furthermore, we perform our duality 
analysis about a fermionic vacuum in which there are two 
and zero d fermions (for an average of one) on the two 
sublattices respectively (the PSG ensures that this choice 
of vacuum actually breaks no symmetries). Alternatively 
stated in more physical terms, the previous analysesiiii^ 
were carried out about a background of fermions occupy- 
ing a half-filled band of single-particle states, while our 
duality analysis assumes that the electrons are primarily 
in paired singlet valence bond states. The consistency 
of our approach with the modulations observed in STM 
studies of the vortex lattice''*' (discussed above) , and the 
absence (so far) of the staggered current patterns near 
vortices predicted by the SU(2) gauge theory^S niay ten- 
tatively be regarded as experimental evidence in support 
of the approach advocated in I and the present paper. 

While our duality analysis is best controlled when the 
fermionic S = 1/2 quasiparticle excitations of the super- 
conductor are fully gapped, the structure of our theory 
does allow these quasiparticles to acquire a gapless nodal 
structure in a superconductor, without a strong quali- 
tative impact on any of the phases of our theory. Upon 
decreasing 5 from deep in the superconducting state, such 
nodal quasiparticles can annihilate in pairs across a re- 
duced Brillouin zone boundary in a superconductor at 
S = Sf, as has been described in previous workiSi^'^^-^o. 
Alternatively, one can imagine that the nodal quasipar- 
ticles persist in the superconducting phase right down 
to the Mott insulator-to-superconductor transition: the 
theory of I will clearly not apply to such a transition, 
and obtaining the correct theory for this case remains an 
important open problem. 
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